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Abstract 

A newly-proposed parity-violating nucleon-nucleon interaction based on effective field theory is 
studied in this work. It is found that at low energy, i.e., Slab < 90MeV for l S Ji Po transitions 
and £i a b < 40MeV for 3 Si~ Pi and 3 Si- 3 Pi transitions, the parity- violating observables can be 
completely specified by a minimal set of six parameters. It contains five low-energy constants, 
which are equivalent to the Danilov parameters, and an additional parameter that characterizes 
the long-range one-pion exchange and is proportional to the parity-violating pion-nucleon coupling 
constant h\. Selected observables in two-nucleons systems are analyzed in this framework with 
their dependence on these parameters being expressed in a nearly model-independent way. 
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I. INTRODUCTION 



Study of the parity-violating (PV) nucleon-rmcleon (NN) interaction, V pv , and its asso- 
ciated nuclear PV phenomena begins with the report by Tanner [l] shortly after the parity 
violation was confirmed experimentally in 1957. The first clear evidence is found a decade 
after by observing a non-zero circular polarization in the 7-decay of 181 Ta |2|. Although 
quite a few PV observables have been measured later on in various nuclear systems, ranging 
from simple two-body scattering to heavy nuclear reaction, our current understanding of nu- 



5 



for reviews of this 



clear parity violation is still far from complete (see, e.g., Refs. 
field). The major difficulty is two-fold: not only these experiments require high precision 
to discern the much smaller PV signals, but also in theory, the non-perturbative character 
of the quark-gluon dynamics makes a "first-principle" formulation of V PV as yet impossible. 
Despite the difficulties and that the underlying theory, the SU(2) <g> U(l) gauge theory, is 
well-established, the study of V pv is still valuable for two main reasons. First, it is the 
only viable venue to observe the neutral weak current interaction between two quarks at low 
energy. Second, it supplies more information about the nucleon-nucleon (NN) dynamics in 
addition to existing scattering data. 

The phenomenological development of V pv proceeds in a similar fashion as what has been 
achieved in the parity-conserving (PC) NN interaction — starting out from pure phase-shift 
analyses, then parameterizations in meson exchange models, and finally to rigorous effective 
field theory (EFT) formulation nowadays. It is first pointed out by Danilov (3, 0, E3| that, 
at low energy, V PV can be characterized by five S-P scattering amplitudes: Af p , A" p , and 
ynn f 1 S Q - 3 P transitions, At for 3 S'i- 1 Pi transition, and p t for 3 5*i- 3 Pi transition. This 
idea is ^aHsed by Deques and M issimer Q Q t „ an effect.ve version — through 
the Bethe-Goldstone equation — which applies to many-body systems. On the other hand, 
formulations of V pv in terms of meson exchange models can be dated back to the seminal 
works by Blin-Stoyle (3, Q| and Barton The specific form involving one 7T-, p-, and u- 
exchanges, Vqme? then becomes the standard in this field after Desplanques, Donoghue, and 
Holstein (DDH) give their prediction of the six PV meson-nucleon coupling constants, /z^'s 
(m denotes the type of meson and % the isospin), based on a quark model calculation 1 As 



1 As the conventional nomenclature, the "DDH" potential, could be somewhat misleading, it is referred as 
the PV one-meson-exchange (OME) potential, Vqmei in this work. We thank B. Desplanques for this 
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explained in Ref. |3j, Vqme has a close connection to the S-P amplitude formulation, Vg V P , 
at low energy: The amplitude pt contains a long-range one-pion-exchange contribution, and 
the other amplitudes, including the short-range part of p t , are all related to the vector-meson 
exchanges. 

Most of the existing PV observables have been analyzed in the one-meson-exchange 
(OME) framework. However, a consistent constraint of the PV coupling constants is not 
realized yet. There are several reasons. On the experimental side, many data have large 
errors so are not very constrictive; also, these observables in terms of h l m 's are not inde- 
pendent enough to allow a simultaneous determination of these six parameters. On the 
theoretical side, several precise data involve many-body systems; the reliabilities of these 
calculations are questionable. As one can see from, e.g., Refs. (j, lla, llifl , a two-dimensional 
constraint on the particular linear combinations of the isoscalar and isovector couplings 
already shows some discrepancy. Besides these possibilities, one might also wonder if the 
analysis framework, i.e., Vqme; cou ld be the culprit. 

To address the last question, Zhu et al. (3| recently reformulate V PV in the EFT frame- 
work to the order of Q (Q is the momentum scale) . This new framework comes with two 
incarnations: one with pions fully integrated out, Vj w , and the other with dynamical pions, 
V£p T . The pionless version, thought to be suitable for low-energy cases with E^b ^ 100 MeV, 
only contains the short-range (SR) interaction, V^g^, and it is specified by ten low-energy 
constants (LECs). However, as Zhu et al. argued semi-quantitatively, only five of them are 
truly independent at low energy, and can be mapped to the five S-P amplitudes. In the 
pionful version, dynamical pions generate two explicit terms: the leading-order, long-range 
(LR) interaction, V^ LR , due to the one-pion exchange (OPE), and the subleading-order, 
medium-range (MR) interaction, due to the two-pion exchange (TPE). While the 

OPE part is also familiar in Vq^ e , the TPE part has never been systematically treated be- 
fore. By considering vertex corrections to the OPE term, the original formulation by Zhu et 
al. contains an extra next-to-next-to-leading-order interaction, V^ R , whose operator struc- 
ture is thought to be different from others already being specified. This term has recently 
been shown as redundant |3| , so will be ignored in our discussion. Overall, in addition to 
the 10 LECs in the SR interaction, the pionful theory introduces, to 0(Q), two more pa- 



clarification. 



3 



rameters: one with the interaction (/t*) and the other with the pion-exchange current (c w ). 
An important point to note is that, although Vf$ R takes the same form in both V^p T and 
Vj v , the LECs in these two EFT frameworks have different meanings: all the pion physics 
is included in LECs for the pionless version; but it is singled out in the pionful version. 

With the advance of experimental techniques showing promise of PV measurements in 
few-body systems — where reliable theoretical calculations are available — an extensive search 
program to re-analyze PV observables is sketched in Ref. The key of this re-analysis is to 
use the "hybrid" EFT framework, which combines the state-of-the-art wave functions (from 
phenomenological potential-model calculations) and the most general form of y PV (from EFT 
techniques). The immediate goal is to find out whether a more consistent picture of nuclear 
parity violation can be reached among few-body systems. The long-term goal of including 
other precise measurements in many-body systems certainly relies on the previous success. 

This paper takes the first step dealing with two-nucleon systems at low energy. The 
aim is to express the observables, both existing and potentially possible, in terms the EFT 
parameters, and to serve as a part of the database which the complete search program calls 
for. The general formalism is introduced in Sec. El The connection between V PV (both 
Vj v and V^p T ) and the S-P amplitudes is studied in Sec. lIIIl To facilitate a both realistic 
and economic search program, special attention is on the quantitative determination of a 
minimal set of parameters and its applicable energy range. The observables of two-nucleon 
systems are discussed subsequently in Sec. IIV( and a summary follows in Sec. El 



II. FORMALISM 



A fully consistent study of nuclear PV phenomena in the EFT framework requires treating 
PC and PV interactions order by order on the same footing. On the other hand, the "hybrid" 
approach, which combines the state-of-the-art wave functions from phenomenology and the 
general operator structure from EFT, is shown to have quite some success. In this work, we 
follow the latter approach as outlined in Ref. 
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A. Parity- Conserving Potential and Wave Functions 



In the hybrid EFT framework, the unperturbed scattering and deuteron (the binding en- 
ergy E-p = 2.22 MeV) wave functions, \ip) and \ip) v , are obtained by solving the Lippmann- 
Schwinger and Schrodinger equations, respectively, 

(H + V pc Tie)m {±) = E\i;)M , (1) 
(H + V pc ) W) v = -E v \^) v , (2) 

with a chosen high-quality phenomenological potential as V FC . In this work, we use Argonne 
vis (AV18) [20] model exclusively. The model dependence of PV observables on strong 
potentials has been extensively studied in Refs. (21I, 22^ : for most cases, no strong deviation 
from AVI 8 is found. 

Since the PV interaction is small, we treat it as a first-order perturbation. The PV 
scattering amplitude, M, is obtained by the first-order distorted-wave Born approximation 

M=(-)(^|\/ PV |^) (+) . (3) 

~(±) ~ 

The parity admixtures of the scattering and deuteron states, \ip) and \i>) v , are obtained 
by solving the inhomogeneous differential equations 



(E-H -V PC M) (±) = \/ pv |^) (±) , (4) 
{Ev + H + V PC M) V = -V™ \^) v , (5) 

respectively, where the product of the PV potential and the unperturbed wave function serves 
as the source term. We refer more technical details regarding the partial wave expansion, 
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23j], but only mention a subtle point 



phase shifts, and numerical procedures to Refs. 
about the phase convention: We adopt, exclusively, the Condon-Shortley phase convention; 
it is different from the Biedenharn-Rose phase convention which contains an additional phase 
i L for the partial wave of orbital angular momentum L. 

B. Parity- Violating Interaction in Pionless EFT 

In the pionless EFT, the PV interaction is entirely short-range and takes the following 
form in the coordinate space [3] 



PV / 
,SRV 



{ [Ct + (C 2 + C 4 ) r z + + C 3 r. + C 5 r zz ] <r ■ y m+ (r) 



+ 



Ci + (C 2 + C4)r; + C 3 r. + C 5 r 2 



fx • Vm-\r) 



+ (C a - C 4 ) o-+ • ?/ m+ (r) + C 6 cr + • y m _(r) 



(6) 



where A x is the scale of chiral symmetry breaking and related to the pion decay constant 
F n by A x = 4nF n « 1.161 GeV; r. = n • r 2 , r| = (rf ± r|)/2, = z (n x t 2 )*/2, and 
r 22 = (3 t{ r 2 z — Ti ■ r 2 ) / (2 y 7 ^) are the isospin operators; 2 cr ± = <Ji±a 2 and cr x = z cr! x cr 2 
are the spin operators. The spatial operator y m -(+)(r) is defined as the (anti-) commutator 
of — i V with the mass 2 -weighted Yukawa function / m (r) 



y m ±[.r) 



-i V, / m (r)] ± 



— z V , m 2 



47rr 



(7) 



± 



When m — > oo, / m (r) — > 5(r)/r 2 ; the potential thus takes a four-fermion contact form as 
expected. 

While using a Yukawa functional form for f m (r) leads to a similar spatial behavior as the 
conventional Vqme) other choices — as long as they are realized in the context of EFT — are 
also possible. For instance, taking into account the monopole form factors at both the strong 



and weak vertices with a cutoff A m , one obtains a modified Yukawa function 



fm(r) 



m 

4nr 



-A m r 



i + 5 (i 



m \ A 



(8) 



iv, nrp. » oo limit, f m { r ) recovers the "bare" Yukawa form f m (r). We note that in 
I21I a recent and extensive OME analyses of two-body nuclear PV, the authors 



At the A r 
Refs. 

adopt f m {r) instead of the conventional choice f m (r). 

As C 2 and C4 appear as a linear combination C2 + C4 in Eq. (@), Vf$ R contains 11 — 1 = 10 
LECs. After resolving the isospin and spin matrix elements of all allowed two-nucleon states, 
the PV observables depend on the following ten linear combinations of C's and C's: 



The operators t zz and t z a we adopt are different from Ref. Q|. Therefore, the LECs C5, C5, and Cq in our 
definition are greater than their counterparts in Ref. Q| by factors of —2 V6, —2 V6, and 2, respectively. 
Also note that the notation of C'e in Ref. is changed into C§ in this paper, because it is associated 
with a y m - type operator like other C's. 
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• pp : Dpp = d + C 3 + C 2 + C 4 + Cs/v^ and Df = C l + C 3 + [C 2 + C A ] + C 5 /y/6, 

• nn : L>™ = d + C 3 - C 2 - C 4 + Cg/v^ and 5™ = Ci + C 3 - [<5 2 + C 4 ] + C 5 /VE, 

• np | Ti=T/=1 : D»* = d + C 3 - 2 Cg/v^ and 5^ = Ci + C 3 - 2 £5/^6, 

• |t 1= t>=o : A« = Ci - 3 C 3 and £>„ = Ci - 3 C 3 , 

• np \t^t s ■ D w = C 2 - C 4 and D w = C 6 . 

In fact, Vig R is tantamount to the p- and cu-sectors of Vqme ^ one assumes i) m p = m w = m 
and ii) the following relations between C"s and C's: 

= = 1 + Xoj , (9a) 

c 3 c 4 c 5 

7r = 7r = 7r = 1 + Xp, (9b) 

o 3 o 4 O5 

where Xu> an d Xp are the isoscalar and isovector strong tensor couplings, respectively. The 
remaining 11—5 = 6 independent LECs in EFT then have a one-to-one mapping to the PV 
heavy-meson-nucleon coupling constants as 

(Ci,<7 2 )--| (10a) 

z, iu N ni w 

(C 3 , C 4 , C 5 , C 6 ) - -| (/£ fcj, fcj, /£) ? (iob) 

where g x denotes the strong x-meson-nucleon coupling constant. Note that in analyses based 
on Vq^ e , the hp' part is usually ignored because it has the same operator structure as the 
LR OPE interaction, i.e., the h\ part, but a very small predicted value for hp . One might 
be tempted to adopt a different Yukawa mass m = m n for the Cq term so that V^g R bears 
even more similarity to Vqj^ e . However, this is undesirable from the EFT point of view 
because two quite different length scales, l/m p and l/m n , both show up in the "short-range" 
interaction. 
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C. Parity- Violating Interaction in Pionful EFT 

When pions are treated explicitly, the EFT PV interaction, as formulated in Ref. 
contains three parts 3 

V^(r) = ^(r) + V^ R (r) + ^ v R (r) . (11) 
The leading term V™ LR is the familiar PV OPE potential 

^-i V lrW = A CI r* <x + ■ y^{r) , (12) 

x 

with 

6 2 V2 F n ml 2 V2 m N ml ' 
where qa = 1-27 is the nucleon axial vector coupling constant and m, = 139.57 MeV; and 
the second equality follows from the Goldberger-Trieman relation. 
The subleading MR interaction is due to TPE and has the form 4 

x 

+ d?ri<r + . [(1-1/(3^)) yL(r)-l/Syg(r)}}, (14) 



with 



C? = -4V2ng A hl, (15) 

C 2 e * = 3V2ng A hl. (16) 

The Yukawa-like radial functions /^(r) and f£i(r), for generating y^ir) and y^ir) via 
Eq. (j2J), are obtained from the Fourier transforms of 



= v^«r±Z ln / y^3±Z±M ) , (17) 

\q\ y 2m ff J 

H{q) = 2 2 ; 2 ^)> ( 18 ) 



3 As mentioned in the introduction, we ignore the higher-order LR term in Ref. [lfij . since it is shown 
to be redundant 0|. 

4 Some typographical errors in Eq. (121) of Ref. have been fixed in order to produce Eqs. CSlEJi see 
also Ref. |2J|. We thank B. Desplanques and Zhu et al. for pointing this out. 



8 



1 1 


i i 1 i i i i 

A=1.31GeV 




A=1.50 GeV 




A=1.72 GeV 




A=10 GeV — 




A— > oo 


1 / 1 


i i 1 i i i i 


0.5 1 



0.2 
0.1 


-0.1 
-0.2 



1 1 1 1 1 


1 1 1 1 1 1 1 1 
A=1.31 GeV — 




A=1.50 GeV 




A=1.72 GeV - 


-V- 


A=10 GeV 




A— > oo 




\ -1.58 r 




oc e 


1 1 1 1 1 


i 



r (fm) 




0.1 



r (fm) 



10 



\ 1 


1 1 1 1 1 1 1 1 
A=1.31GeV 




A=1.50 GeV — 




A=1.72 GeV 




A=10 GeV - 




A— > oo 




\ -1.47 r 




\ oc e 


1 1 1 1 1 


i — i — 



1 



r (fm) 



10 



Figure 1: The r- weighted Yukawa-like radial functions r /^(r) and r f^(r) in the two-pion- 
exchange potential. 

respectively. In Fig. HI the plots of r f^ir) and r f£^(r) are shown with several dipole 
cutoff factors (A 2 — 4m 2 ) 2 / (A 2 + g 2 ) 2 — including the bare case, i.e., A — > oo — introduced in 
the Fourier transforms. As one clearly sees, the short distance behaviors are quite cutoff- 
sensitive, especially for /^(r) since L(q) diverges logarithmically as \nq/m n . On the other 
hand, the long-range tails, roughly decrease like e~ L58r and e~ L47r , track well with e - 2m ^ r = 

e -1.41r 

The way we define Cg , Cf^ and is handy for the bookkeeping purpose; this gives 
V-i,lr an d ^i ? mr the same formal structure as the corresponding parts — as hinted by the 
subscripts — in Vf^ R (also, all the Yukawa functions f m (r), f n (r), /^(r), and/^(r) have the 
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same limit when m, m n — > oo). However, this does not imply these pion PV constants are 
comparable in magnitude with LECs C's and C's. Comparing Eqs. (|T3| fT5| fT6|) . one sees 
C|g smaller than Cg roughly by an order of magnitude. By Eqs. (|10a,[ llObj) and assuming 
all the 7T-, p-, and cu-coupling constants approximately the same, one estimates C's smaller 
than Cg roughly by a factor of m 2 p /m\ ~ 30; for C's, due to the tensor couplings, Eqs. (jHal 
l9b|) . the suppression can be less. If the above assumptions are not too far off, we can roughly 
conclude that C|g and the LECs, C's and C's, are of the same order, and all of them smaller 
than Cg by an order of magnitude. This observation is consistent with the power counting 
scheme that both the TPE and SR terms are of the same higher order than the OPE one. 
But, more definitive answer should still be sought from experiments. 

We also emphasize that the TPE term V^ R , Eq. (JHj) , being used here is not the full PV 
TPE potential. According to Ref. Q|, it only contains the singular part of the TPE, and 
the rest regular terms are still effectively included in the short-range interaction. In this 
sense, V^r m f ac t depends on the chosen regularization scheme. However, as long as the 
most general operator structure is maintained and one does not try to fit the data over a 
large energy range, a consistent analysis should be scheme-independent — this will become 
clear in the next section. 

D. Setup and Parameters 

In the following sections, various PV observables in two-nucleon systems will be ana- 
lyzed by V^p T . While the full results are realized in the pionful EFT framework, the ones 
correspond to the pionless EFT framework can be easily read by simply retaining only the 
part from Vf^ R — with the notion that the LECs in this case effectively include all the pion 
contributions. Although the choice of the Yukawa mass parameter in Vi$ R is arbitrary, we 
use the p meson mass, m = m p = 771.1 MeV, since it has an easy connection to the meson- 
exchange picture. For the convenience of presentation, these calculations will be referred as 
the "bare" calculations, because all Yukawa functions in V^p T are not modified by any form 
factor. 

The above results will be checked against existing calculations in the Vqme framework. 
This is done by applying the relations Eqs. (|Ha| l9b | llObl I10a,[ andHBJ) to V^p T and ignoring 
all the TPE contribution; from now on, we call this procedure OME-mapping (OME-m). 
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Table I: Sets of strong parameters used for OME-mapping. 



g-K 




9 P 


9u> 




K p 




DDH-best [3j 13.45 
DDH-adj. [22] 13.22 




2.79 
3.25 


8.37 
15.85 




3.70 
6.10 


-0.12 
0.0 


Table II: Sets of weak parameters 
systems, h° p + h l p + h 2 p /VG = -22.3 


used for OME-mapping (in units of 10" 
instead of —24.8 as shown in this table, 


~ 7 ). Note the 
will be used 


it for pp 
21|. 


hi 


h° 
n p 


hi 


t 


h° 


hi 


K 


DDH-best [15J 4.56 
DDH-adj. [21, 22] 4.56 


-11.4 
-16.4 


-0.19 
-2.77 


-9.50 
-13.7 


-1.90 
+3.23 


-1.14 
+1.94 


0.00 
0.00 



For numerical estimates, the strong parameters are taken from Ref. jb| and the weak ones 
are set to be the "best-guess" values of DDH [15]. This set is labeled as "DDH-best" in 
Tabs. □ and HU 

In order to _ with Re,, Q, as pointed ont previous,, one has to use the 
monopole-modified Yukawa functions instead. For this matter, we perform a parallel set of 
calculations using f m (r) with the cutoff parameters chosen to be A p = A w = 1.50 GeV and 
A,r = A 27r = 1.72 GeV. These calculations will be referred as the "mod" calculations. While 
using f p ,u,(r) in V^ R does not contradict the EFT framework, using f^{r) and f^i.f) does 
not seem fully consistent with EFT. For the OPE part, this is less a problem because form 
factors only suppress the LR OPE slightly. On the other hand, the validity of adding form 
factors to the TPE part in EFT needs further justification. Thus, our calculation for this 
TPE part should only be understood as showing a qualitative feature in cases where form 
factors are built in. For numerical results in this set of calculations, the strong parameters 
are taken from the Bonn model [3] and the weak ones are the fitted results of Refs. j^l], I3- 
This set is labeled as "DDH-adj." in Tabs. HI and iHl 

III. S-P AMPLITUDES AND DANILOV PARAMETERS 

A PV potential with 11 (10 LECs plus hi) undetermined parameters certainly poses a 
formidable challenge — how can we gather sufficient data and do reliable theoretical analyses 



11 



of them? A substantial reduction of the LECs is proposed in Refs. [13] by building the 
connection between the EFT framework and the S-P amplitude analysis, which is pioneered 
by by Dando, fl Q B, and .tended by Deques and M issimer .atee on QUR. 
The main idea of this reduction goes like the following. 

For low-energy PV phenomena in which only the S-P mixings contribute substan- 
tially, the observables can be expressed by five independent PV scattering amplitudes: 
v pp ,nn,np(} So ^ 3 Pq), u( 3 Si ^ Pi), and w( 3 Si ^ 3 Pi). From the last section, we know 
each amplitude due to the Vfg R part is a linear combination of the corresponding D and D 
with the coefficients being determined by the matrix elements of y m+ and y m - , respectively. 
An important observation comes from that the matrix elements (y m +) and (y m -) are equal 
in the zero range approximation (ZRA) with the absence of the NN repulsive hard core. 
This causes D and D always appear in a (D + D) combination and work actually like one 
energy-independent LEC. Therefore, the number of LECs can be reduced to 5 which cor- 
responds to the number of independent S-P amplitudes. While this sounds an attractive 
idea, however, neither ZRA nor the absence of NN hard core are physically realized. In 
order to put this reduction scheme of LECs along with the whole search program proposed 
in Ref. [Tj| on a firmer ground, we try to answer a series of questions which have not been 
addressed: 

1. For realistic cases, this 10-to-5 reduction can still work as long as the matrix elements 
of y m+ and y m - have (almost) the same energy dependence, i.e., the condition 

^^- = R{E)^R, (19) 

{Urn—/ 

is satisfied. Therefore, we try to determine R and its range of constancy over E for 
each S-P amplitude. 

2. When higher partial waves become important, the S-P analyses are no longer valid. 
Thus, it is necessary to determine the maximum energy for this proposed search pro- 
gram to work. 

3. In Ref. Q, it is also proposed that for E lah < 100 MeV, the pionless EFT should 
work. However, this requires the contributions from V_i LR and caTl be effectively 
included in V^ R . We try to justify whether this condition can be met. 
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4. At the end of this section, we determine the zero-energy S-P amplitudes, the so-called 
Danilov parameters, in terms of the PV parameters in V^P^. They will be the actual 
parameters used to express the PV observables in the next section. 

According to the definitions by Desplanques and Missimer [lfl|. the S-P amplitudes are 
calculated by the following formulae: 



m N ( 3 P 


\V py 1 
1 ^eftI 






ip ( 3 P 


(T ■ f 


I'So) 


r(i + z77)| |r(2 + z?7)| 



v~ D v + v + D v + v 27T Cf , (20a) 



u 



w 



-m N 


( 1 -Pi|^eftI 3 'S'i) 


ip 




u~ D u 


+ u + D u , 


-m N 


(^il^rl^i) 


ip 


^P 1 \ ( T + -f\^S 1 ) 



(20b) 



= w-D w +w + D w + w*C%+w 2 "C% r , (20c) 

where all the amplitudes are functions of energy, p is the two-nucleon relative momentum, 
and the factors are chosen to reproduce the normalization and limiting behaviors of Refs. 
^ 1^1 . Note that for the v amplitude, an extra factor, which is 1 when the Sommerfeld 
number 77 = 0, is introduced in order to completely subtract the Coulomb effect for pp 
scattering at threshold. 

Since the total cross section is proportional to the imaginary part of the forward scattering 
amplitude, we concentrate on Im[f,w, w]. 

The top panels of Fig.[2]show the energy dependences of x~ [x G v pp ^ nn ^ npi u, w) type S-P 
amplitudes, proportional to (y m _), up to E\ ah = 100 MeV. Because the form factors suppress 
the short-distance contributions, the results from the "mod" calculations are consistently 
smaller than the "bare" ones. Another noticeable feature is the plots of v pp , v nn , and v np 
overlap, and the tiny difference is mostly due to the small charge-dependent interaction built 
in AV18. 

The bottom panels of Fig. [2] show the percentage deviations of x + /x~ from their zero- 
energy values x + (0)/x~(0), i.e., 

X ,X ~ z+(0)/*-(0) ' [2L) 
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Figure 2: Top panels: the energy dependences of the x~ type S-P amplitudes. Bottom panels: the 
percentage deviations of x + jx~ from their zero-energy values x + (0)/x~(0), where the dotted lines 
mark the 10% level. 

In case the energy dependences of x + and x~ keep the same, Ax + /x~ remains zero. There- 
fore, this quantity is a measure of the departure from the perfect scaling, Eq. ()19|) . which a 
strict 10-to-5 reduction scheme requires. As these plots show, the deviations all grow with 
energy in the positive direction. For the v and w amplitudes, the scaling actually works very 
well — though not perfect — up to = 100 MeV where the deviations are still less than 10% 
for the "bare" case. For the u amplitude, however, the 10% tolerance for scaling deviation 
can only allow go as high as 40 MeV. It might come as a surprise why the u and w 
amplitudes have such different behaviors, since both involve the same 3 Si wave and the 1 P\ 
wave (for u) does not differ from the 3 Pi one (for w) dramatically. The answer is due to the 
tensor force, by which a true distorted 3 Si-wave acquires some -D-wave component. From 
the ratios 



(*P 1 \tr--r\*Si,L = 2}/( l P 1 \(T_-r\ 3 S 1 ,L = 0) = -V2 , 



( 3 Pi|<x + • f| 3 Si, L = 2)/( 3 P 1 |ct + ■ f| 3 Si, L = 0) = 1 



(22) 
(23) 
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Figure 3: The percentage corrections induced by the D P transitions to the S-P ones for the x 
type amplitudes, where the dotted lines mark the 10% level. 



one learns that the D-wave component is more enhanced in the u than the w amplitude, 
and this is the main cause of the difference. A mock calculation by pretending the ratio in 
Eq. J22J to be l/y/2 as in Eq. verifies that the scaling of u + /u would then be similar 
to w + /w~. The "mod" calculations generally show larger deviations, so the ranges within 
which the scalings work have to be reduced. This can be understood from that the difference 
of y m + and y m - is a term involving the gradient on the wave function. As the form factors 
suppress the short-range contribution, the longer-range part of the wave function which has 
a larger gradient thus gets a bigger weight; this leads to an enhancement of the deviation. 

The importance of the DP transitions is illustrated in Fig.OJ where their percentage 
corrections to the S-P transitions are shown for the x~ type amplitudes. Given a 10% 
tolerance for corrections, one can conclude that the S-P transitions dominate up to = 
90 MeV for the v type amplitudes (it can be higher for the pp case because of the Coulomb 
repulsion) and up to = 40 MeV for both the u and the w types (note that the w 
type contains two contributions from 3 D 1 - 3 P 1 and 3 D 2 - 3 P2). The "mod" calculations do 
not change these conclusions too much. It should be stressed that these D P amplitudes 
have quite different energy dependences compared to the S-P ones. Therefore, although 
the scaling behaviors found above for the S-P amplitudes can apply to higher energies, the 
10-to-5 reduction should be limited by the prerequisite of the S-P dominance. 

In Fig.Hl whether the OPE and TPE contributions can be effectively integrated out 
and lead to a purely pionless short-range potential is studied by examining the scaling 
behaviors of one- and two-pion amplitudes x 7T,27T with respect to x~ . As one can see, all 
the deviations increase with energy in the negative direction. In the "bare" calculations, the 
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Figure 4: The percentage deviations of x 7T,27T jx~ from their zero-energy values af^(0)/x~(0), 
where the dotted lines mark the 10% level. 

TPE amplitudes scale very well with the SR one. Allowing a 10% deviation, the scaling 
works as Ei ah reaches up to more than 100 MeV for the v amplitudes and 60 MeV for the w 
one; both limits are above the S'-P-dominant region. On the other hand, the scaling works 
quite poorly — only up to £i a b = 20 MeV — in the "mod" calculations. This can be easily seen 
from Fig.|H where the modified two-pion Yukawa-like functions differs substantially from the 
bare ones at short distances. Since we have already mentioned the inconsistency of adding 
form factors to the TPE potential in such an ad hoc fashion in EFT, the curves labeled by 
"mod" for the TPE part should not be taken too seriously. 

The most noteworthy information in Fig. SI comes from the observation that it is almost 
impossible for the OPE amplitude to scale with the SR one, even within a modest energy 
range, say 10 MeV or so. From threshold to 40 MeV, in which the S-P dominance holds 
for the w amplitude, the scaling deviation increases to 50%. Thus, this re-confirms the old 
wisdom that it takes two parameters — one for the SR and the other for the LR terms — to 

n 

characterize the physics of the 3 Si~ 3 Pi transition p|. In this sense, the applicability of the 
pionless EFT to nuclear PV is extremely limited, only within a very narrow energy range 
near threshold. 

One might be temped to think that, as this bad scaling is due to the huge difference 
between the Yukawa mass scales: m SR = m p in V^ R and m L R = m n in V^ LR , then, 
choosing a smaller mass scale in V^ R should make the pionless EFT work. The most 
straightforward choice is msR = — as long as the energy being considered is much smaller 
than m n — by which the w~ and amplitudes become identical. In this Fig.na 
shows (the "mod" calculations are very close the the "bare" ones since the form factors do 
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Figure 5: The similar analysis as Fig.0 for the x~ and x + amplitudes with the Yukawa mass 
parameter in V"_- SR chosen to be m = m^. The dotted line marks the 10% level and it is not 
visible for the central bottom panel. 

not cutoff the LR interaction too much), the scaling between x + and x~ becomes extremely 
bad, and this completely spoils the nice reduction scheme of 10-to-5 LECs. For the v and 
w amplitudes, the 10% scaling deviation only allows E\ a \> < 5 MeV; for the u amplitude, it 
is even worse (the divergence at E\ ah ~ 50-60 MeV is because Im[w~] crosses over the zero 
value). The reason is again that, by making the effective range of V^ R longer, the difference 
between (y m +) and (y m -) gets enhanced as the long-range part of the wave function, having 
a larger gradient, gets a bigger weight. Therefore, the price to pay for making the pionless 
EFT work is to keep all the 10 LECs as independent parameters. Clearly, the more economic 
choice would be keeping the OPE part explicitly and having the 10-to-5 reduction work by 
choosing ttj-sr = m p for V^ R . 

To summarize the discussions so far, we conclude that at low energy, i.e., Ei ah < 90 MeV 
for 1 S'o- 3 Po transitions and Ey ah < 40 MeV for 3 Si- 1 Pi and 3 5i- 3 Pi transitions, an PV po- 
tential formulated in EFT requires at least 6 parameters: 5 LECs plus C§ (x. h\. One may 
think that this conclusion makes V^p T equivalent to Vq^ e , which also has six PV parame- 
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ters: 5 heavy-meson- nucleoli couplings (ignoring h^) plus h\. However, this statement is not 
true in general, because the six-parameter EFT framework only works under the assump- 
tion of the S-P dominance, but the OME framework is not limited by this requirement. 
Furthermore, the OME framework implies some prescribed relationships, Eqs. (fTTal l9b| HOal 
and HObj) . between C's and C"s in the EFT framework; unfortunately, these relationships 
can not be justified without going to high energy. 

At last, we come to the determination of the Danilov parameters, Af p '" n,np , A*, and p t , 
which will serve as the 5 independent LECs. The relations between Danilov parameters and 



the zero-energy S-P amplitudes are given by 

v NN ' (0) « -af N ' e l «S£ (°» A™' , (24a) 

u(0) « -a t e* {6lpl {0)+63S1 m X t , (24b) 

w S r(0) « -a* e^iW+Wo)) p t , (24c) 

where a denotes the corresponding scattering length and 5(0) the zero-energy phase shift (in- 
cluding the Coulomb contribution) |lfl| j. The notation "wsr" means the OPE-subtracted 3 S\- 

3 Pi amplitude. Using the values af = -7.8064 fm, a™ = -18.487 fm, < p = -23.7318 fm, 
and a t = 5.4192 fm, we get the dimensionless Danilov parameters 

m N \PP = 5.507 x 10~ 3 (5f + 0.789 - 1.655 C^) , (25a) 

m N \™ = 5.796 x 10~ 3 (5™ n + 0.792 D™ n - 1.648 <5 2 27r ) , (25b) 

m H X* = 5.778 x 10~ 3 (IX* + 0.809 D" p ) , (25c) 

m N \ t = -1.462 x 10~ 3 (D u - 2.230 D w ) , (25d) 

miV p t = 3.108 x 10~ 3 (D w + 0.604 D w - 1.771 C7f ) ; (25e) 

for the "bare" case; and 

m N Af = 3.268 x 10~ 3 (5f + 0.849 D p v p - 1 .260 C 2 2,r ) , (26a) 

m N \™ = 3.809 x 10" 3 (D™ + 0.853 D™" - 1.237 <5 2 2 ") , (26b) 

m N X^P = 3.772 x 10" 3 (D^ p + 0.871 D™ p ) , (26c) 

mjv A t = -0.867 x 10~ 3 (D u - 2.425 D u ) , (26d) 

m^pt = 2.003 x 10~ 3 (D w + 0.664 D w - 1.586 C$*) , (26e) 

for the "mod" case. 
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In these expressions, one sees that the scaling factors between D's and D's are not 1 
as anticipated by the ZRA with the absence of hard core. The most striking case, which 
amounts to 2.230 (or 2.425 for the "mod" case), happens for the 1 Pi- 3 Si transition. This 
in fact shows the attractiveness of the hybrid EFT treatment for nuclear PV: one gets the 
strong dynamics right without having to go to a higher order in Q, where the proliferation 
of needed PV parameters is totally undesirable. 



IV. PARITY- VIOLATING OBSERVABLES IN TWO-NUCLEON SYSTEMS 

Having determined a minimal set of PV parameters, i.e., m N \ p v ,nn,np , m N \ t , m N p t and 
Cg, to describe the low-energy PV phenomena, we will use them in this section to express 
the PV observables which have been or will be measured in two-body systems. As analyses 
of these observables have been quite extensively discussed in the Vqme framework, we refer 
most of the details which also apply to the EFT framework to literature and only highlight 
the new results and the comparison with the old framework. 



A. A 1 ^ in pp scattering 

The "nuclear" total asymmetry for pp scattering 
Af(E) 
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(27) 



is defined through the Coulomb-subtracted, forward (9 = 0) scattering amplitude A4(E,Q), 
where the subscript pair S'M' S , SMs denotes the final and initial two-body spin states, 
respectively. As the Coulomb scattering amplitude, M c , diverges at the forward angle, 
the total asymmetry (with full 4n angular coverage) can only be well-defined after this 
singular piece is subtracted: M. = M — M c . One should note that A p £ is not a quantity 
an experiment directly measures; a theoretical correction is needed to fold an experimental 
result into A^(E) (see, e.g., Refs. [^[ill] for more discussions). 

Currently, there are two low-energy data points at 13.6 MeV and 45 MeV which give 
Af = -(0.93 ± 0.21) x 10- 7 Q, Q| and -(1.57 ± 0.23) x 10" 7 Q, respectively. These 
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Table III: Analysis of A pp decomposed in partial waves. Each entry denotes the multiplica- 
tive coefficient for the corresponding PV coupling constant. The full result is the sum of every 
"entry x coupling" in the same row. The last column "OME-m" shows the numerical value of A p £ 
by performing the OME-mapping to the EFT result with the PC and PV parameters specified in 
Tabs. HI and El 





^o-^o (xlO- 3 ) 


1 D 2 - 3 P 2 (xlO' 3 ) 


1 L> 2 - 3 F 2 (xlO- 3 ) 


OME-m 

(xlO- 7 ) 


d pp 


d pp 




D pp 


yfp 




d pp 


yfp 




13.6 


bare 


-1.980 


-2.476 


4.010 


-0.005 


0.006 


-0.012 


0.001 


-0.001 


0.002 


-0.971 


mod 


-1.398 


-1.617 


1.885 


-0.004 


0.004 


-0.010 


0.000 


-0.001 


0.001 


-0.960 


45 


bare 


-3.686 


-4.476 


7.026 


-0.122 


0.132 


-0.243 


0.027 


-0.033 


0.064 


-1.746 


mod 


-2.582 


-2.868 


2.842 


-0.089 


0.094 


-0.180 


0.019 


-0.023 


0.050 


-1.662 


221 


bare 


-0.069 


-0.073 


0.112 


-2.749 


2.618 


-3.784 


0.164 


-0.388 


0.633 


0.426 


mod 


-0.046 


-0.043 


0.015 


-1.888 


1.678 


-1.636 


0.086 


-0.270 


0.420 


0.853 



supersec 



e the earlier, less accurate experiments at 15 and 45 MeV which yield — (1.7±0.8) x 



one measurement at 221 MeV, yielding +(0.84 ± 0.29) x 10~ 7 |3iJ [3j|; it is motivated by 



10- 7 (aJ, lasj and -(2.31 ± 0.89) x lO" 7 |2S [sj, respectively. At hig her energy, there is 



y. At hig her e 

7 0, H; ^ if 



the theoretical prediction that this would uniquely determine the PV p-exchange coupling 
constant hf = h° p + h 1 p + h 2 p /V& in V^ E Q- 

The EFT analysis of A P l for the low-energy data points is tabulated in Tab. IIII1 Appar- 
ently, the observables are dominated by the S-P transition. Using the Danilov parameters 
obtained in the last section, we find that 

Af (13.6 MeV) = -0.449 m N Xf + (-0.035 - 0.088 x 10" 3 (bare), 

or - 0.445 m^Af + (-0.032 D pp - 0.157(7^) x 10~ 3 (mod); (28) 

Af(45MeV) = -0.795 m^Af + (-0.329 - 0.395 C 2 ?) x 10~ 3 (bare), 

or - 0.771 m^Af + (-0.276 D pp - 0.813 Cf) x 10~ 3 (mod). (29) 

The correction terms, enclosed in parentheses, are in general quite small except for the 
TPE parts in the "mod" calculations, so we conclude that these two data only depend 
on one single parameter, m N \ pp . In fact, the theoretical prediction for the ratio of 
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Af (45MeV)/Af (13.6 MeV) « -0.795/ - 0.449 = 1.77 (or -0.771/ - 0.445 = 1.73 for the 
"mod" case) agrees very well with the experimental result ~ —1.57/ — 0.93 = 1.69 (discarding 
errors). Furthermore, one can see from the comparison between the "bare" and "mod" results 
that even though m N \ p s p and m N \ pp are defined in different models, the expressions for A pp 
in terms of them are almost model-independent. This justifies the advantage of using the 
Danilov parameters instead of the LECs C's, C's, etc. 

Even though it is doubtful that the V^p T of order 0(Q) is sufficient for analyzing the 
high-energy datum at 221 MeV, it is nonetheless interesting to see how the analysis turns 
out to be. The result, also shown in Tab. IIII[ indicates a quite different feature. The 
D-P amplitude becomes the most dominant contribution, with the D-F one also being 
non-negligible. Both amplitudes have their own scaling factors between D pp , D pp , and Cf^ 
components different from the S-P amplitude. While it is not clear if the D P and D-F 
amplitudes can be completely specified by D pp , D pp , and C% n , they certainly can not be 
uniquely fixed by the only high-energy measurement. 

By OME-mapping the EFT results in Tab. EH our calculations are checked with litera- 
ture. The "bare+DDH-best" results are consistent with works such as Refs. 



3. |22,Q,Q,Q, 



given that different strong potential models are used. The "mod+DDH-adj." results agree 



well with Refs. 
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, the small difference is because we do not use a different mass and a 
cutoff factor for the u meson. 

It is worth to point out that the analysis by Carlson et al. j^l, which is based on a OME 
framework with two independent parameters h pp and h^f, claims a good fit to both low- and 
high-energy data. Unfortunately, due to the lack of more high energy data, it is currently 
impossible to verify this fit along with its important dynamical assumptions — the monopole 
form factors and big isovector tensor coupling Xp — within the EFT framework. On the other 
hand, their fitted PV uNN coupling constant, h pp = /i° + h^, is only marginally consistent 
with most hadronic perditions and needs further clarification. For these issues, we refer to 

n 

Ref. (40j for more discussions. 

Finally, we turn our attention to the TPE contributions, which have not received extensive 
study and are left out while we do the OME-mapping in Tab. Mil By writing out C% n in 
term of h\ explicitly and assuming the DDH best value for h\, one sees the asymmetry in 
the "bare" case increased by ~ 70% for the 13.6 and 45 MeV data points and ~ 60% for the 
221 MeV one. In the "mod" case, the increases are ~ 30% and ~ 20% for low- and high- 
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energy cases, respectively. Although these corrections seem quite big, one should remember 
they are not the full TPE corrections, as part of the TPE contributions are buried in the SR 
interaction. But since and have the same power counting, it is not unnatural to 

expect individual terms of similar magnitude. Another remark concerns the general trends 
that the TPE enhances the asymmetry and it is the the low-energy cases that get more 
boost than the high-energy ones. These points have also been noticed in Ref. 4£j| , where 
part of the TPE contribution is accounted for by formulating it as a p-resonance. 



B. cj)n P and Pn P in neutron transmission through para-hydrogen 



It was first pointed by Michel |41| and later on by Stodolsky [42, 43] that nuclear parity 
violation can be studied through low-energy neutron transmission, where the whole process 
acts like optics. The observables could be a spin rotation, n , about the longitudinal axis 
(assumed to be z) for the transversely-polarized neutron, or a net longitudinal polarization, 
P n , that an unpolarized neutron beam picks up when traversing through medium — the latter 
is equivalent to the asymmetry in cross section for the longitudinally-polarized neutron 
scattering, A 1 ^. These quantities per unit length (assuming the target is uniform), d(p n /dz 
and dP n /dz, can be related to the PV forward scattering amplitude M(E, 0) by 

^ = -^N Re(M +z (E, 0) - M_,(E, 0)) , (30) 
^ = -^N lm(M +z (E, 0) - M_ Z (E, 0)) , (31) 

where k is the magnitude of the neutron momentum, N is the target number density, and 
the subscript ±z denotes the direction of neutron polarization. 

For a thermal neutron beam, E n ps 0.025 eV, transmitting through liquid par a- hydrogen, 
N = 0.24 x 10 23 /cm 3 , the EFT analysis of d<p n /dz and dP n /dz is tabulated in Tab. ED At 
thermal energy, the magnitude of d P™ p / dz is about four orders of magnitude smaller than 
d(p^ p /dz. When neutron energy is further decreased, d(j)™ p /dz stays constant, but dP™ p /dz 



or low-energy 
about the 



drops as y/E n ; therefore, the spin-rotation measurement is more feasible 
neutrons. This trend is consistent with the argument made by Stodolsky 
elastic scattering. It is also pointed out in Ref. (43j that exothermic processes, i.e., inelastic 
exit channels, can possibly result in a non- vanishing dP n /dz at zero energy. However, it 
is not the case for np scattering, where the only exothermic reaction, np — > dj (will be 
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Table V: Analysis of (I) d4hf{th..)/dz in rad/m and (II) d Pn P (th.)/d z in 10 4 /m decomposed in 
partial waves. See Tab. II I II for the explanation of tabularization. 





1 5 - 3 Po(xlO- 2 ) 


3 Si- 1 Pi (xl0~ 2 ) 


3 Si- 3 Pi (xlO" 2 ) 


OME-m 

(xHT 7 ) 


D7 




D u 


D u 


D w 


D w 


ci 


°6 


(I) 


bare 


1.169 


I AAA 


-0.186 


0.083 


0.214 


0.355 


0.286 


-0.628 


6.711 


mod 


0.821 


0.943 


-0.120 


0.049 


0.152 


0.229 


0.284 


-0.363 


5.149 


(II) 


bare 


4.805 


5.939 


0.175 


-0.079 


-0.202 


-0.334 


-0.269 


0.591 


2.165 


mod 


3.378 


3.878 


0.113 


-0.047 


-0.143 


-0.215 


-0.267 


0.341 


-2.220 



discussed in Sec. IIVD|1 . does not lead to a total asymmetry, as remarked in Ref. |4j 

In this case, all three different S-P amplitudes come into play, and the results in terms 
of Danilov parameters and Cq are 

d<j>f (th.) 



d z 



m/rad 



0.286 CI + 2.500 m n \ n J> - 0.571 m N \ t + 1.412 m N p t + (0.000) (bare); 

or 0.284 Cq + 2.500 m N \™ p - 0.571mA + 1.412 m N p t + (0.000) (mod). 

(32) 

Because this process is close to zero energy, it is not a surprise that the Danilov parameters 
work extremely well (almost no error). Also, this expression is model-independent, no matter 
it is for the "bare" or the "mod" calculation. 

By OME-mapping the EFT results, the "bare+DDH-best" value, d(j^{t\\.)/dz ~ 6.71 x 
10 -7 rad/m, is about 20% smaller in magnitude than an early prediction using the Paris 
potential and with a different sign. Thus, we confirm the assertion of Ref. about 
the sign problem in Ref. As for the "mod+DDH-adj" value, ~ 5.15 x 10 -7 rad/m, it 
agrees well with Ref. If these numerical estimates are not too far off, the plan of doing 
such an experiment aiming at a 2.7 x 10~ 7 rad/m precision at the Spallation Neutron 
Source (SNS) will certainly provide a valuable data point. 

The TPE contribution enters through the 3 S'i- 3 Pi transition. Because Cg and C% n have 
the same sign, Tab.lVlsuggests that the TPE reduces the OPE contribution which dominates 
the above OME-m estimates. The correction is about —15% for the "bare" calculation, and 
—10% for the "mod" calculation. This ~ 10% correction is consistent with the qualitative 
power-counting argument that the TPE contribution is smaller than the leading OPE one 
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Table VI: Analysis of P™ p (th.) decomposed in partial waves. See Tab. II I II for the explanation of 
tabularization. 





l S - 3 P mix. (xl(T 3 ) 


V- 1 P 1 mix. (xl0~ 3 ) 


OME-m 

(xlO- 7 ) 


d7 


D np 


D u 


Du 


bare 


-0.751 


-0.935 


2.166 


-0.980 


0.247 


mono 


-0.525 


-0.607 


1.391 


-0.580 


0.520 



by an order of magnitude. 



in np — ► c?7 and A\ in 7c? — > np 

Low-energy radiative neutron capture mainly involves the lowest-order electromagnetic 
transitions. For np — > dj, it is Ml for the PC part, and El for the PV part. Since the 
total cross section is dominated by the 1 So-wave scattering, the non-zero circular polarization 
takes an approximate simple form as 





\ 3 Po) 


+ V 




\Ei\ 


1%) 




1%) 



(33) 



where the double bar "||" denotes the reduced matrix element. In this case, the observable 
depends on the 1 Sq- 3 Po and the deuteron T)- X P\ admixtures. It is important to recognize 
that we rely on the Siegert theorem j^, through which the E\ operator is related to the 
charge dipole operator C\, to calculate the El matrix elements. This manipulation not only 
implicitly includes most 0(v/c) meson exchange currents, but also imposes a AS = spin 
selection rule as shown in Eq. (|33|) . 

For thermal neutron, the EFT analysis is tabulated in Tab. IVII Although the observ- 
able P™ p is not determined directly by the scattering amplitudes — instead, by the parity 
admixtures which do not have a trivial relation to the scattering amplitudes in general — the 
Danilov parameters still do a good job 

P"P(th.) = -0.161 m JV A" p + 0.670 m N \ t + (0.005 D n u p + 0.019 D u ) x 10~ 3 (bare); 

or - 0.161 m N Xl P + 0.669 m N \ t + (0.004 D^ p + 0.016 D u ) x 10" 3 (mod). (34) 

The reason is mainly because the process occurs at a very low energy and the deuteron is a 
loosely bound state; both are not far from the zero-energy limit. 
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The OME-mapping gives the "bare+DDH-best" value P 7 p (th.) = 2.5 x 10 8 which agrees 
well with a recent calculation j^] and is also consistent with pre-80's predictions, e.g., 
Refs. Q, Q, Q Oh around (2-5) x 10" 8 , (see Ref. Q for a summary). For the 
"mod+DDH-adj." value P™ p (th.) = 0.52 x 10~ 7 , we have an excellent check with Ref. Q- 

Historically, the first — of P 7 (.h.) done by the Leningrad group reports a 
result —(1.3 ±0.45) x 10 [54], which not only exceeds most theoretical predictions by two 
orders of magnitude but also has an opposite sign. The follow-up experiment does correct 
the sign problem; however, the published result P" p (th.) = (1.8± 1.8) x 10 -7 [55] still has too 
large an error bar. In order to circumvent the difficulty of measuring a circular polarization, 
the inverse process, the asymmetry Ajf - in deuteron photo-disintegration, j d — > np, can 
be a good alternative. By detailed balancing, = P" p if all kinematics are exactly 
reversed. One can show that, for photon with energy of 1.32 keV above the threshold, 
A| d (1.32keV+) = P 7 np (th.). 

As demonstrated in several theoretical works (3,0,0,0,01, the asymmetry A]f gets 
larger when approaching the threshold, but on the other hand, the total cross section gets 
smaller. There are two data points reported in 80s: (2.7 ± 2.8) x 10 -6 at P 7 = 4.1 MeV and 
(7.7 ± 5.3) x 10" 6 at P 7 = 3.2 MeV Q, 0]; though they qualitatively justify the statement 
above, the precisions are too low to be of use. With various groups showing interests of new 
measurements (e.g., Ref. 0]), ^ * s important to decide the best photon energy (should not 
be too far from the threshold for larger asymmetry) to be employed. 

An important point to note for this particular observable is that, unlike the case for 
neutron transmission, the T>- X P\ admixture has an important contribution so that the model 
dependence is worrisome. The situation is most clear when c omp aring with other semi- and 
non-local potential model calculations. As shown in Refs. (3,01, the CD-Bonn and Bonn-B 
calculations give predictions two times larger, and the Bonn calculation even enhances by an 
order of magnitude. The difference is mostly due to the large variations of these models in 

n 

the 1 P\ channel at short distances [48]. In this sense, a well-determined X t , if ever possible, 
can be used in a reversed way to constrain strong potential models. 
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Table VII: Analysis of ^4" p (th.) decomposed in partial waves. See Tab. II I II for the explanation of 
tabularization. 





3 Si- 3 Pi mix. (xl(T 3 ) 


V- 3 P 1 mix. (xl0~ 3 ) 


OME-m 

(xl(T 7 ) 


D w 


D w 


ci 


°6 


D w 


D w 


ci 


°6 


bare 


-0.108 


-0.185 


-0.133 


0.321 


-0.066 


-0.103 


-0.139 


0.193 


-0.506 


mod 


-0.076 


-0.118 


-0.132 


0.172 


-0.048 


-0.069 


-0.138 


0.128 


-0.486 



D. A-y in np — » 

By the same approximation as in the above subsection, the photon asymmetry in ftp 
c?7, A™ p , which is defined through da±{6) / 'dfl oc 1 ± A™ p cos 9, 5 can be expressed as 



A np 



■y/2 



(V\ 


\Ei\ 




+ V 


( 3 ^i 


\El\ 


I 3 S1> 


(V\\Mi\ 


1%) 



(35) 



Unlike the case for P" p , it is the 3 Si- 3 Pi and T>- Z P\ admixtures that contribute in this case. 

For thermal neutron, the EFT analysis is tabulated in Tab. IVII1 Again, the Danilov 
parameter m N p t pretty much summarizes the short-distance physics 

A" p (th.) = -0.272 CI - 0.093 + (0.003 C% w ) x 10~ 3 (bare); 

or - 0.270 & % - 0.093 m N p t + (0.004 C^) x 10~ 3 (mod). (36) 

-5.06 x 10~ 8 for the "bare+DDH-best" case and 



The OME-mapping values, A" p (th.) = 
— 4.85i< IO7 8 for the "bare+DDH-adi." case, are consistent with existing predictions, e.g. 

for the former and Refs. 



for the latter, 



Refs. 3, Q, Q, B, 0, 0, Q 

respectively. 

The TPE contributions change the above OME-m results somewhat. Their corrections 
to the OPE contributions are —13% and —8% for the "bare" and "mod" cases, respectively. 
This is similar to the neutron spin rotation case, and consistent with a recent calculation [24]. 

The great interest of measuring A™ p is mainly because it is dominated by the OPE in 
the Vqme framework. This can also be observed from the above EFT analysis: If one 
assumes the natural size of m N p t /CQ ~ 0.1, the OPE contribution then dominates the SR 



5 From this expression, one clearly sees J dftda + (9) = J dQd<j-(9). This confirms the earlier statement 
that A 1 ^ vanishes at zero energy, even if an exothermic process exists. 
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one by a factor of 30 or so. 6 One of the outstanding puzzles in nuclear PV is the difficulty 
of accommodating the extremely small upper limit on h\, set by the 18 F results j^j], with 
hadronic predictions and other nuclear PV experiments. The NPDGamma experiment 
currently running at the Los Alamos Neutron Science Center (LANSCE) and will be at SNS 
later on, aims to reach an ultimate sensitivity of 5 x 10~ 9 . Results from this experiment will 
certainly improve the long-existing value: (0.6 ± 2.1) x 10 and hopefully resolve 

the h\ puzzle. 

Concluding this section, we shall make an important remark about the PV meson ex- 
change current (MEC) effects which manifest in electromagnetic processes such as radiative 
neutron capture being discussed here and in Sec. IIV CI Although the Siegert theorem 
alleviates much of the problem regarding the calculations of transverse electric multipole 
operators E/s due to MECs, there is no easy simplification when the transverse magnetic 
multipole operators Mj's are concerned. Furthermore, the Siegert theorem only applies to 
MECs which are constrained by the continuity equation; for other transverse MECs, their 
effects to E/s have to be added separately. 

In Ref. there is indeed such a transverse MEC, which can not be accounted for by 
gauging the PV interaction, and it introduces a new PV constant designated as c n . This 
MEC takes the following form in the configuration space 

je n (x]Xi } x 2 ) =-i — jp% ( r i+ r 2-) [o-2-fo-i x Va.<y (3) (aj-aJi)] 

p-mjr r 

x r (l + m w r) + (1^2), (37) 

where r± = r x ±ir y , and r = \x\ — x 2 \. Compared with the dominant part of the PV OPE 
MEC, the so-called pair current, 

j^ivix; x u x 2 ) = - - * - (n ■ t 2 - t{ tI) [<ti 5 {3) {x - Xi)] 
2y2m N 

+ (1^2), (38) 



Airr 



the matrix element (j 57r ) roughly scales with (jVpair) by a factor (— i V x /A x ) = k/A x , assum- 
ing (r) ~ l/m n for typical nuclei. Therefore, for the radiative processes considered in this 



6 We stress that this argument is based on naturalness. Without further experimental confirmation, one 
should still keep other possibilities open. 
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work, where the photon energy k is just a few MeV so that k/A x < 1%, the contribution of 
jgjt is negligible. Hence we do not have to include this extra PV constant cv in the current 
search program at low energy. 

V. SUMMARY 

In this work, we study the newly-proposed search program for nuclear parity violation 
based on the effective field theory framework [18]. It is found that, in oder to completely 
describe the nuclear PV phenomena at low energy, a minimal set of six parameters is needed. 
By low energy, it means Ey ah < 90 MeV for processes involving 1 Sq- 3 Po transitions and 
E\ab ^ 40 MeV for ones involving 3 S\- 1 Pi and 3 S\- 3 Pi transitions. The six parameters to be 
determined phenomenologically are the five dimensionless Danilov parameters: m N X pp,nn ' np ', 
m N X t and m N p t , and the long-range one-pion-exchange parameter Cg, which is proportional 
to the parity- violating pion-nucleon coupling constant h\. 

The two-body parity- violating observables being studied in this work are summarized as 
following: 

Af (13.6 MeV) w -0.45 m N X p s p , (39) 
Af (45 MeV) « -0.78 m N X pp , (40) 
^-07(th.)| rad/m « 0.30 Q + 2.50m JV Af - 0.57 m N X t + \Alm N p t , (41) 

P" p (th.) « -0A6m N X^ p + 0.67 m N X t w (1.32 keV+) , (42) 
Af(th.) w -0.27 (7^ - 0.093 . (43) 

Because A^ p (13.6 MeV) and A^ p (45MeV) essentially determine the same quantity, m N X pp , 
these equations only serve as four constraints — if precise data can all be obtained for the 
three listed neutron experiments. In order to have at least two more linearly-independent 
equations, other experimental possibilities have to be explored. In few-body systems, where 
reliable theoretical analyses can be performed, the candidate reactions include pd, nd, pa, 
n a etc. — just to name a few. Currently, there are a published datum for p a: A^ Q (45 MeV) = 
— (3.3 ± 0.9) x 10~ 7 [70], and an ongoing experiment of thermal neutron spi n ration in liquid 
helium, 0^ Q (t/i.), at the National Institute of Standard and Technology |4a| . In this respect, 
existing calculations of these PV five-body processes should be updated. Because alpha 
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particle is a tightly bound state such that nucleons inside have larger momenta, whether the 
S-P dominance — the cornerstone of this six-parameter analysis — can still hold should be 
carefully examined. On the other hand, low-energy reactions involving d or t might suffer 
less the problem. However, in order to motivate new experiments, updated theoretical works 
are indispensable. 
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